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We present a significant improvement over our previous calculations of the cosmic string contribu- 
tion to cosmic microwave background (CMB) power spectra, with particular focus on sub-WMAP 
angular scales. These smaller scales are relevant for the now-operational Planck satellite and addi- 
tional sub-orbital CMB projects that have even finer resolutions. We employ larger Abelian Higgs 
string simulations than before and we additionally model and extrapolate the statistical measures 
from our simulations to smaller length scales. We then use an efficient means of including the 
extrapolations into our Einstein-Boltzmann calculations in order to yield accurate results over the 
multipole range 2 < £ < 4000. Our results suggest that power-law behaviour cuts in for £ > 3000 
in the case of the temperature power spectrum, which then allows cautious extrapolation to even 
smaller scales. We find that a string contribution to the temperature power spectrum making up 
10% of power at ^ = 10 would be larger than the Silk-damped primary adiabatic contribution for 
£ > 3500. Astrophysical contributions such as the Sunyaev-Zeldovich effect also become important 
at these scales and will reduce the sensitivity to strings, but these are potentially distinguishable by 
their frequency-dependence. 



I. INTRODUCTION 



Observations of the cosmic microwave background 
(CMB) radiation and the large-scale distribution of 
galaxies indicate that cosmic structure was seeded in the 
very early stages of the universe [ll y], consistent with 
the inflationary paradigm. However, the datasets leave 
room for significant effects due to the presence of cos- 
mic strin gs a t subsequent times [3|-|6|. These strings (see 
Refs. |7h10|| for reviews) are particularly important in 
that they are predicted by many physically motivated 
inflation models, including brane inflation models in the 
context of string theory [ll|, [l2| but also models rooted 
in Grand Unified Theory (GUT) [H. Such models also 
predict other types of cosmic defect, including textures 
and semilocal strings [14|-[l7|. and their CMB signals are 
[18l - l22| also of great interest. 



In previous works we calculated cosmic string CMB 
temperature and polarization power spectra using field 
theory simulations of the Abelian Higgs model [23, [2J| 
and used the results to fit models with both string- and 
inflation-induced anisotropics to CMB data [3- It was 
found that the data (which was primarily the WMAP 
third- year release [23]) favoured a model with a fractional 
string contribution to the temperature power spectrum 
at multipole moment £ = 10 of /lo = 0.09 ± 0.05 with 
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the spectral tilt of the inflation-induced primordial per- 
turbations being n^ — 1.00 ±0.03. The latter was in con- 
trast to the inflation-only result of rig = 0.95lto'oi9 01 
and showed, along with Refs. [J, [23], that the inclusion 
of topological defects could readily allow Ug > 1 under 
those data, thanks to a parameter degeneracy found in 
Ref. [2^. 

In the present article we present a significant improve- 
ment in the CMB power spectrum predictions from string 
simulations, with particular emphasis on small angular 
scales. While the numerical simulations of the Abelian 
Higgs model in Ref. [23| yielded CMB results cover- 
ing multipoles £ Ki 2 —^ 1000, the full range of angu- 
lar scales relevant for WMAP data, they did not have 
the dynamic range to accurately investigate finer angular 
scales. Higher multipoles are now of great interest thanks 
to good coverage by sub-orbital CMB experiments [271- 
l3l| , plus the imminent arrival of full-sky Planck data [32| 
{£ < 2500), and we have therefore focused on yielding ac- 
curate string predictions for £ < 4000. 

Our ability to now include smaller scales stems par- 
tially from improvements in computational facilities and 
the carefully-chosen initial conditions that we now em- 
ploy in our Abelian Higgs simulations. But more impor- 
tantly, we have made developments in the method used 
to yield the CMB power spectra from the statistical mea- 
sures of the energy-momentum distribution in our simu- 
lations: the unequal-time correlation functions (UETCs) 
[isL l33l |. We now model and extrapolate the measured 
UETCs to smaller scales in order to mimic results from 
larger simulations and have developed an efficient means 
of including UETC results at these extrapolated scales 
into our CMB calculations. 



The need for the above extrapolation can be under- 
stood as fohows. The width of a cosmic string is micro- 
scopic (perhaps at the GUT scale) while their separation 
at times of cosmological importance is of order the Hub- 
ble distance and hence it is not possible to simultaneously 
resolve both scales in numerical simulations during the 
epochs when strings would have impacted on the CMB. 
We can solve this problem by using scaling 34]: when 
the horizon is greater than about one hundred times the 
string width, we observe that the strings enter a late-time 
attractor solution in which statistical measures of the 
string distribution, when measured in horizon units, are 
constant in time. For example the average string length 
in a horizon volume divided by the horizon size is ~ 50. 
Similarly, the UETC measurements are independent of 
time once scaled by the horizon size and hence scaling 
provides knowledge of their values at times of cosmolog- 
ical importance. However, the measured UETC power is 
attenuated on scales close to the string width and scal- 
ing is broken for such length scales. Hence the resultant 
CMB results would show too little power at very small 
angular scales: the sharp changes in temperature caused 
by the Gott-Kaiser-Stebbins (GKS) effect [H^ would 
have been smeared out by having been effectively sourced 
by strings whose width was around one- hundredth of the 
horizon size in the post-recombination era. For our pre- 
vious work this was of little actual concern since those 
articles focused upon WMAP scales {£ < 1000). It is, 
however, necessary to extrapolate the UETC results to 
sub-string scales in order to obtain accurate results from 
field theoretic simulations for much higher multipoles. 

Other approaches to calculating CMB perturbations 
from strings include using simulations of the Nambu- 
Goto (NG) type l38l440l and the unconnected segment 
model (USM) (il. l4l| - |43| . Nambu-Goto simulations can 
yield a greater dynamic range than field theoretic simu- 
lations but must still invoke scaling to yield CMB results 
over a wide range of scales. Furthermore the network cor- 
relation length used in the initial conditions appears to 
persist and may limit the reliable resolution of Nambu- 
Goto simulations [4J,|43]. The USM represents the string 
network as a stochastic set of moving "sticks" , which dis- 
appear at an appropriate rate in order to give the chosen 
string scaling density. While it has no true dynamical 
content, it is computationally cheap and offers the fiex- 
ibility to choose the coarse-grained network properties 
[4,42, 46 48]. 

The advantage of the Abelian Higgs model is that it 
includes the small-scale physics near the string width, 
which has a non-negligible impact on the string dynam- 
ics [4^, [5^: energy from the strings is converted into 
massive gauge and Higgs radiation. In NG simulations 
this decay channel is not included and the string length 
density is significantly higher, with the long strings be- 
ing converted into small loops that would then decay via 
gravitational radiation (although this process is not ac- 
tually simulated). With the extra decay channel in the 
field theoretic simulations, decay via gravitational radia- 



tion would be less important, and this fact significantly 
changes the predictions for gravitational wave observa- 
tions. In the case of the CMB on the other hand, it 
is the long strings that are important, and the key dif- 
ference between the simulation results is the inter-string 
separation. A potential disadvantage of a field theory 
simulation is that computational constraints require the 
string width to be artificially increased in order to keep 
it above the simulation resolution, but we carefully show 
that this does not significantly affect the UETCs and 
therefore the CMB power spectra results. 

In Sections |lT] and IIIII we detail our methods, includ- 
ing an overview of the UETC approach, our field theory 
simulations, the tests of scaling that we employ and the 
sub-string extrapolation. We exhibit the resulting CMB 
power spectra in Section ITVl and give conclusions in Sec- 
tion El 



II. METHOD OVERVIEW 

As already noted, the basis for our string simulations 
and CMB calculations is our previous work: Ref. |23| . 
referred to hereafter as BHKU. We refer the reader to 
that article for the full details of our approach but we 
present the essential information in what immediately 
follows and highlight the improvements made for the cur- 
rent article in the next section. 



A. UETC approach 

In CMB power spectra calculations for inflationary 
models with cosmic strings, the inflationary contribution 
is essentially uncorrelated with the cosmic string con- 
tribution. This is because the inflationary perturbations 
are laid down by an independent field 60 e-foldings before 
the strings are formed, and because the complex string 
dynamics rapidly destroys correlations with earlier times 
(see Sec. IIIIC[) . As a result we may write the total spec- 
trum i{£ + l)Ce as a sum of two independent spectra and 
take the cross-correlation term to be negligible: 



C/^inf I /ns 
I — Lf + Co 



(1) 



where C™^ is the infiationary contribution and Cf''^ is the 
string contribution. We can use standard methods [51[ 
to determine C"^ and it is therefore upon Cf"'' that this 
article is focused. 

Physically the strings cause CMB anisotropics by cre- 
ating perturbations in the space-time metric, which are 
roughly of the same order as Gfi, where /i is the string 
mass per unit length, G is the gravitational constant, 
and Gjj, < 10~^ to be consistent with current observa- 
tions [3|. These inhomogeneities in the metric then lead 
to perturbations in the matter and radiation that them- 
selves evolve and infiuence the strings, but we can neglect 
this back-reaction since the resulting perturbations of the 



strings would then result in changes only of order (Gn)'^ 
to the metric. 

In order to determine the string contribution to a two- 
point correlation function, such as the CMB temperature 
power spectrum, we are required to solve a set of lin- 
ear differential equations of the following form, in which 
the string energy-momentum tensor components act as 
source terms Sa- 



Vacik,a,p,...)Xa{'k,T) = S'c(k, r). 



(2) 



Here Xa is the quantity of interest and Xa its Fourier 
transform, while Vac is a differential operator, depen- 
dent upon the cosmic scale factor a, the background 
matter density p and similar quantities. Our notation is 
such that T is the conformal time and k is the comoving 
wavevector. The homogeneous version of this equation 
{Sc = 0), which corresponds to the inflationary case, can 
be solved by standard codes and therefore, in principle, 
we may use a Green's function Gacik,TQ,T) to give the 
power spectrum at conformal time tq for the string case 
via: 



'XaiK To)lb*(k, To)) ^J Jdrdr' Gacik, T„.T)gUk, Tq, t) 

'5,(k,T)5:(k,r')).(3) 



Hence the data required to calculate such two-point cor- 
relation functions are the two-point unequal-time corre- 
lators of the string energy-momentum tensor components 
MM: 



i/a6(fc,T,T') = (5a(k,T)^,*(k,T'; 



(4) 



Note that statistical isotropy implies that Uab is not de- 
pendent on the direction of k and further that it is real- 
valued. 

Significant simplification in the form of Uab niay be 
made using the scaling property, briefly mentioned in the 
introduction. Under scaling, any statistical measure of 
the spatial distribution of strings scales with the horizon 
size, which is just t in comoving coordinates. For exam- 
ple the comoving length-density of string is a/r^ , where a 
is a dimensionless constant. The existence of this attrac- 
tor was predicted by Kibble 1341 and has been confirmed 
in Nambu-Goto simulations [3^, |4^, [S^] (with some as- 
sumptions about the decay of loops) and in Abelian Higgs 
simulations [23, 49, 53]. We will present further evidence 
in support of scaling in our results section. 

The power of scaling is that it enables us to write Uab 
in terms of a function of just two variables: 




\tIt'). (5) 



;7afc(k, t,t') = 



Here 0o sets the energy-scale of the problem and con- 
verts between the scaling spatial distribution of string 
and the distribution of energy. For example, in the case 



of the Abelian Higgs model (see next section) it is the 
vacuum expectation value of the scalar field. The co- 
moving simulation volume V appears here because we 
define our Fourier transform so that it leaves dimensions 
unchanged: 



l(k) = -^ / d^x X(x) e 



— zk-: 



(6) 



The UETC scaling function Cab can be seen to allow data 
to be taken for very large k at small r and r' and then 
used to provide information about small k at large r and 
t'. This is critical for the present calculations, as noted 
in the introduction. 

Further power of the UETC scaling functions derives 
from their functional form: they decay for large and small 
time ratios t/t' , and for large k^frr' (small scales), while 
causality constrains their form at low k^rr' . This cru- 
cially means that we need study a scaling network for 
only a relatively short range of times and in a limited 
simulation volume. However as noted in the introduc- 
tion, small length scales become more important when 
considering small angular scales (see next section). 

While our approach is to measure Cab from cosmic 
string simulations, and therefore find Uabi we do not in 
fact then use Eq. ([3]). Instead we decompose Cab into a 
sum over products as [SJ]: 



Cab{ky 



• r/r') 






\nCna{kT)Cnb{kT'). (7) 



The problem then breaks down to: 



X,(k,r)X*(k,r')) = '^Y.^r,^a{k.r)ir{Kr'), 



where: 



V 



i:{k,To)= / dTgab{k,To,T) 



Cnbikr) 



(8) 



(9) 



In practice we do not calculate this integral via a Green's 
function, but instead apply a modified version of CM- 
BEASY [55[ to determine the CMB power spectrum con- 
tribution from the coherent active source Cnbj \f^- 

In principle there are 55 possible UETCs between the 
4 scalar, 4 vector and 2 tensor degrees of freedom in the 
energy-momentum tensor T^j/, but thanks to statistical 
isotropy and energy conservation, we in fact need only 
to measure 5 scaling functions [53. The scalar UETCs 
that we calculate involve projections from T^^ that, via 
Einstein's equations, directly source the two Bardeen po- 
tentials I57I: 



cS 
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T„ 



oKfjiKfi-^ 7 



m'^n-^ mn- 



(10) 

(11) 



From these projections there are 3 independent UETC 
scahng functions that we must measure: 

5|(k,r)5|*(k,r')) = -^^Cfi(fcV77^,r/r'), (12) 

(^|(k, r) ^r (k, r')) ^^^Cf.ikVTT^, t/t'), (13) 



5|(k,r)5|*(k,T')W^^C|2(fc^/77,r/r'). (14) 



Then in the tensor case we project out the two tensor 
degrees of freedom (see BHKU), which we denote as S'^^ 
and S'^'^. We then determine: 



S^\k, t) 5Ti'(k, r')) = (^^"(k, r) ^T2*(k, r')) (15) 
= 2- 



/)4 1 ^ 

=2=-CT(fc%/77,r/T') 



where the factor of 2 is present to ensure that C^ matches 
the definition of Ref. [I^ ■ 

Fohowing Ref. [1^, we make a change in the scahng 
function definition when considering the two vector de- 
grees of freedom from Tq,, S^^ and 5^^ (ggg bHKU). 
Specifically we pull out a factor of k'^rr' from the scaling 
function such that: 



5^i(k, t) 5^1 (k, r')) = (5^2(k, r) 5^^ (k, r' 



(16) 



This definition is motivated by the fact that covariant 
energy-momentum conservation requires that the vector 
UETC varies as fe^ at small k while the other UETCs 
that we measure tend to a constant value as fc — > 0; and 
it is desirable for all UETC scaling functions to have the 
same super-horizon properties. 

However, it should be noted that the above discussion 
requires a small change because scaling is broken near the 
time of radiation-matter equality Toq , since Tgq is a second 
dimensional scale which enters the problem. We hence 
are required to take UETC data in both the radiation and 
matter eras — although the matter era data dominates 
the CMB results — and we then use interpolation in or- 
der to model the transition. Scaling is also broken as the 
Universe enters the current accelerating phase, causing 
the string density to decay, which we model by partially 
suppressing the stress source. 

For more details of the solution of the linearized 
Einstein-Boltzmann equations in the presence of sources 
seeRefs. IH,!!!. 



B. Field theoretic simulations 

For computational speed, we simulate local cosmic 
strings by solving the classical field equations of the 
simplest theory that contains them: the Abelian Higgs 



model. In the notation of BHKU this has Lagrangian 
density: 



C 



1 



A 



F,,F^'^ + {D,^r{D^cl^)~-{\<l^\^-4>l) 



(17) 
with D^^ = df_, + iA^ and F^^ = d^A^ - d^A^ while 
e and A are dimensionless coupling constants. In a 
spatially fiat Friedmann- Robertson- Walker (FRW) met- 
ric with scale factor a, this leads to the following Euler- 
Lagrange equations: 



'ci> + 2-j~ D,D,<P = -a^^ (|0|2 _ ^l) 0, (18) 

Foj~d^F,j - -2a^e^Im[cj)*D,c^], (19) 
-diFoi = -2a^e'^Im[(j)*^], (20) 

where the gauge choice Aq — has been made. Here 
we use over-dots to denote differentiation with respect to 
conformal time t, while dt denotes differentiation with 
respect to comoving Cartesian coordinates. 

When simulating this model in an expanding universe, 
the simulations must resolve the comoving string width 
wo/a (i.e. wq is the fixed physical width ~ (J)q ). They 
must also contain at least one horizon volume of co- 
moving diameter 2r. However, these two scales diverge 
rapidly, with their ratio growing as t^ in the matter era, 
while we require ratios > 100 in order for strings to 
scale. Hence with 1024"^ lattice simulations, which are 
the largest that are practical with our available facilities, 
we cannot study a scaling network using these equations 
for long enough to measure the UETCs up to sufficiently 
high t/t' ratios. Furthermore, the increase in this ratio 
with computer time icpu varies as tcpu^ and hence the 
returns from much larger outlays are minimal. In BHKU 
we therefore proceeded by allowing temporal variations 
in the coupling constants A and e: 



A = 



A, 







a2(l-s)' 

eo 



a^ 



such that the comoving string width now varies as 






(21) 
(22) 

(23) 



where s is the string width control parameter. If s = 0, 
then the factors of a on the right-hand-side of Eqs. P^ - 
()20p are removed and the string width remains constant 
in comoving coordinates, while s = 1 gives the normal 
dynamics of the model. Note that the above dependen- 
cies preserve the ratio A/2e^, which we set to be unity for 
the simulations described in this article: the Bogomolnyi 
limit [13. 

Simply evolving the above dynamical equations (Eqs. 
1181 and I19p while varying A and e will not preserve the 
constraint Eq. (|20p . Hence the BHKU method is to vary 
the model action with respect to the fields whilst allowing 



for the temporal dependence in A and e, which then yields 
a consistent set of dynamical and constraint equations: 
Eqs. ([TSl) and ([20)) in addition to a modified form of 
Eq. dUD: 



Fc 



oj -r 2(1 - s)-Foj - d^F.j = -2a^e^ Im[(f)*Dj(l)] . (24) 



However, for s y^ 1 the action is no longer a 4-scalar 
and hence there is a breach of covariant energy conserva- 
tion. Since it is the energy-momentum tensor that seeds 
the cosmological perturbations which result from strings, 
then this is clearly a potential problem. Fortunately in 
BHKU it was established that the effects on the UETCs, 
and therefore the CMB power spectra, are minimal on 
the relevant scales, but we present additional evidence in 
support of this approach in the next section. 
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III. METHOD REFINEMENTS AND 
INTERMEDIATE RESULTS 

As discussed in the introduction, the use of scaling 
to translate simulation results from GUT length-scales 
to cosmological scales means that any effects present in 
our simulations on scales close to the string width are 
erroneously transfered to scales of order one hundredth 
of the horizon size at the decoupling of the CMB, i.e. 
^ 3 kPc rather than ^ 10^"^^ m (which would correspond 
to the string width if G/i ^ 10^^). As already mentioned, 
this smearing out of the energy density would lead to 
reduced power in the string contribution to the CMB 
power spectrum at small angular scales. 

For our previous work, which was intended to be 
used only with WMAP data, this effect was not antic- 
ipated to be significant, since WMAP only probes scales 
larger than about one tenth of the horizon at decoupling 
[i < 1000). Additionally the small-scale data which then 
existed was not precise enough or on small enough scales 
for the likely inaccuracies to be a cause of concern, unless 
the strings completely dominated the CMB on such scales 
— something that the WMAP data was seen to rule out 
Q. While it is true that at times long after recombi- 
nation, the properties of strings on scales much smaller 
than the horizon can influence WMAP scales, their effect 
would have been minor, as we shall demonstrate in this 
article, which now includes their contribution. 

In order to include small-scale UETC power, we model 
the UETCs and extrapolate the trends seen to sub-string 
scales fSec lIIIDJ) . We have also modified our CMB cal- 
culation method in order to rapidly include these extrap- 
olations (Sec. IIIIE[) . However, we additionally employ 
larger simulations than in the past and employ different 
initial conditions, both of which enhance our ability to 
study the scaling epoch and improve our measurements 
of the UETC scaling functions. We hence discuss our 
new results for measures of scaling and of the UETCs 
themselves, before going on to discuss the inclusion of 
small scales. 



FIG. 1: Results for ^ from simulations using our new ini- 
tial conditions compared to those using BHKU initial con- 
ditions, which yielded an offset scaling law. Statistical un- 
certainties determined from 3 realizations are denoted using 
the shaded regions (1 — ct dark, 2 — cr light) and the linear 
best-fits over the period between the short vertical lines are 
indicated by the dashed lines. The results come from 1024'', 
s — 0, Ax = 0.5(/>(7 simulations in the case of the new initial 
conditions and 768^, s = 0.3, Ax = 0.42950^^ for the BHKU 
ones. 



A. Tests of scaling 

As has already been alluded to, for strings to scale 
their width must be much less than their separation, 
which is of order the horizon size. This statement may 
be made more precise by introducing the (comoving) net- 
work length scale ^, defined by: 



e=\^ 



V 



(25) 



where L is the comoving string length^ in the simulation 
volume V. We observe in our simulations that accurate 
scaling sets in at ^ > 40u' whereupon ^ ^ 0.3t. 

However, in BHKU we were limited to an Abelian 
Higgs simulation size of 512'^, which is about the min- 
imum required to study a scaling network of strings. We 
therefore were forced to be content with initial conditions 
that yielded not ^ oc r but ^ oc (r — tj=o), as shown in 



^ We measure L by detecting the lattice grid-squares around which 
the phased has a net winding (using the gauge-invariant method 
of Ref. 16111 ) and then we approximately reconstruct the string 
path as a collection of perpendicular segments of length Ax pass- 
ing through these squares. We then apply the Scherrer-Vilenkin 
correction factor of n/6 |62|| in order to approximately account 
for over-estimate due to representing a smooth path by perpen- 
dicular segments (see Ref. [50| for more discussion). 



Fig. [H were rj=o is negative. That is, at early times 
string decay occurred more quickly than scaling would 
have predicted and f increased rapidly, but then stabi- 
lized to an approximately constant d^/dr once ^ > AOw 
(where w « 0q but is slowly decreasing with r in the 
s — 0.3 simulation for which BHKU-style results are plot- 
ted). As can be seen in the figure, this rapid increase 
aids the creation of an approximately scaling network for 
simulations that are causally limited to small times, al- 
beit scaling with (t — tj=o) rather than r. Fortunately, 
at times of cosmological interest we have r ^ t^=o s-nd 
therefore this offset is not actually relevant. Furthermore, 
when T was replaced by (r — r^^o) in the calculations of 
the UETC scaling functions, the equal-time scaling func- 
tions C{kT, 1) were seen to be approximately constant in 
T, as is required for scaling. 

Despite the success of the (t — r^=o) replacement when 
using the BHKU initial conditions, with the larger 1024^ 
simulations available for the present work we noted a 
slow drift in the d^/dr value, even at late times. Fur- 
ther, there was a slow rise in the magnitude of the 
Tj^o-corrected equal-time scaling functions. This is ex- 
pected at some level since the dynamical equations con- 
tain the quantity: 
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a 



n 

r 



(26) 



where n is unity for the radiation era and 2 for the mat- 
ter era. The ratio of this damping scale r/n to the net- 
work length scale ^ changes with time if t^=o t^ and 
at T ~ \t(=o\ with negative tj=oj the system has a lower 
than asymptotic value of r/n^ — it is being damped too 
heavily. We hence now employ initial conditions that 
yield (effectively) t^=o = 0, although now we do not see 
scaling until long after the causal run-time limit of our 
previous 512'^ simulations. 

In principle, the initial conditions for the fields are 
set by their fluctuations (quantum or thermal) at string 
formation, which have a finite, microscopic, correlation 
length (see e.g. [63] )■ Fortunately the string network is 
observed to relax to scaling for a wide variety of initial 
conditions and so it is not important to model the ini- 
tial fluctuations precisely. In BHKU we employed initial 
conditions that were designed to model a vacuum phase 
transition at the end of inflation: each lattice site was 
given an independent phase for 0, with \(f)\ = (po, and the 
initial time Tstart was set roughly equal to the lattice spac- 
ing Ax. The gauge field and the canonical momentum (jj 
were set to zero. However, that set of initial conditions 
means that the simulations begin with Tstart ^ w and 
this was seen in BHKU to be responsible for the initial 
rise in ^. 

In this work we begin with </) as a Gaussian random 
field with correlation length l^ such that l^ ~ Tstart ^ w. 
We are free to choose the two-point auto-correlation func- 
tion of this Gaussian random field, subject to the con- 
straints that it is zero outside the causal horizon and 
that its Fourier transform (the power spectrum) must be 



AH 
AH 
NG 



Lt^/V 
Lr^/V 



0.255 ±0.018 

15 ±2 

37.8 ±1.7 



0.285 ±0.011 
12.2 ±0.96 
28.4 ±0.9 



TABLE I: Numerical results for the network length scale ^ in 
horizon units and the string length density L/V normalized 
to the horizon size. AH indicates that the results are from 
the present Abelian Higgs simulations while NG indicates the 
results are from the Nambu-Goto simulation of Ref. 44] (see 
Sec. IIV A2[) . which yield string densities approximately 2.5 
(radiation) and 2.3 (matter) times greater. (Note the NG 
results quoted includes only infinite strings, although loops 
do not contribute significantly to our figures !50|.) 



non-negative. For simplicity we relax the first constraint 
slightly and opt for: 



P{y) = ^j'^x' 



'(x-y)=Foexp(V/2Zl). (27) 



This decays so rapidly outside the horizon that it is effec- 
tively causal with horizon size ^ l^. We hence generate 
a field in Fourier space with spectrum P and random 
phases for each k-mode, before transforming it to real 
space to become our initial (j) configuration. This gives 
{(t>*(f)) ~ Pq, which we take as <j>q since the field should be 
close to its vacuum except at lattice sites very near to the 
strings. Further, since Hubble damping at t ^ w is too 
weak to rapidly relax the fields into a network of string, 
we now employ diffusive (first-order) evolution until a 
time Tdiff in order to achieve that. We then have three 
parameters Tstart, Tdiff and l^ which we may vary in or- 
der to achieve our goal of t^=o = 0, although each test 
of a point in this 3-dimensional parameter space is very 
computationally-expensive. 

While in principle this t^=o — >■ optimization should 
be performed for each value of s, the insensitivity of the 
string dynamics to changes in s (as noted in BHKU) im- 
plies that this computationally-costly process can be per- 
formed only once. However, to keep the same settings for 
these 3 parameters at all s values, we must use a lower 
lattice spacing Aa; at higher s in order for strings to still 
be resolved {wo/a'^ > 2Aa;) at the end of the simulation 
T = NAx/2, which is the causal time limit^. This implies 
that the accessible t/t' range narrows as s is increased. 

The results for i^ from these new initial conditions can 
be seen in Fig. [l] for simulations in the matter era using 
s = (with a = 1 , A = 2 and e = 1 during diffusive evolu- 
tion and Ax = 0.5(/)g" ). Further, we show from the same 
simulations examples of the equal-time scaling functions 



^ We continue the simulations slightly beyond the strict causal 
limit, since they do not feel the periodic boundary conditions 
immediately, but the further reduction in string width is not 
sufficient to greatly affect the reliability of the simulations. 





FIG. 2: The equal-time scaling functions C{kT, 1) from 1024^, 
s = simulations with a matter-dominated FRW background 



for 5 times, as indicated in the legend in units of 



'0 



Re- 



sults shown are the average of 3 realizations, with the esti- 
mated statistical uncertainties for t = 296(^^^ indicated by 
the shaded regions. 



FIG. 3; The equal-time scaling functions C{kT, 1) from siiu- 
ulations with a matter-dominated FRW background at three 
values of s. Results shown are the average of 3 realizations, 
with statistical uncertainties for the s = 0.5 case indicated 
by the shaded regions. Results are plotted for r = 150(^,7^1 
which is the start of the period when we take UETC data for 
the s = Q case. 



C{kT, 1) in Fig. [21 which provide a scale-dependent test 
of scaling and is the quantity that is most important for 
the CMB calculations. These results exhibit scaling on 
large and intermediate scales to a good degree of accuracy 
over a ratio in conformal time of 2. A breach of scaling 
is evident at small-scales, where statistical uncertainties 
are highly suppressed, but as noted in the introduction 
this is expected due to the proximity of the associated 
length scales to the string width. We discuss this issue 
and our solution to it in Sec. IIIIDI Numerical results for 
the string length density under scaling are also given in 
Table U (see also Sec. lIVA2|) . 



B. Dependence of the equal-time scaling functions 
upon s 

In addition to ensuring that the simulations are ex- 
hibiting scaling, we must also check that the use of s < 1 
in our dynamical equations does not introduce significant 
systematic errors in our UETC scaling function results. 
This is most readily illustrated using the equal-time case 
again, with results shown in Fig. [3] We can resolve no 
dependence upon s, except for at high fcr, which is again 



due to the proximity of the corresponding length scales 
to the string width. Since the affected scales are the very 
ones for which we apply the extrapolation mentioned in 
the introduction, then this effect is of no concern (see 
Sec. IIIIDJ) . For completion, however, the s-dependence 
on these scales can be explained as follows. The higher 
s, the more rapid the reduction in comoving string width 
during the Hubble phase and therefore the attenuation 
of small scale power (see Sec. IIIID[) manifests itself at 
higher kr. Hence the s — case yields lower results at 
the highest-plotted kr values. 

Since s = simulations enable the greatest range in 
t/t' under scaling, while accurately matching the dy- 
namics seen at higher s values, our final CMB results 
will be based upon s ~ simulations and we will limit 
our remaining discussion in this article to simulations at 
this value of s. 



C. UETC scaling function results and decoherence 

Despite our improved initial conditions, we are only 
able to study the system when it is scaling accurately for 








FIG. 4: The UETC scaling functions C{ky/TT, t/t') in the matter era from 1024^ simulations with s = 0, averaged over 3 
realizations. The raw data is highlighted by the lighter central region, with extrapolations to more extreme t/t' values and to 
lower k\^TT' indicated by the darker regions. The vertical axis indicates \C\ and it should be noted that the cross-correlator is 
negative near its peak. 



conformal time ratios t/t' ~ 2 (corresponding to phys- 
ical time ratios ^ 4 in the radiation era and ~ 8 in the 
matter era). As can be seen in Fig. |4] this is sufficient 
to map out the important region of the UETCs, and to 
permit extrapolation in t/t' as explained below. It can 
be seen that the auto-correlator scaling functions peak 
for T — t' and decay for unequal times, with decay oc- 



curring for t/t' ratios that deviate only slightly from one 
if ky/TT' is large but more slowly on super-horizon scales. 
In the cross-correlation case (7^2 , this is broadly true but 
the peak is noticeably offset on super-horizon scales. 

This behaviour can be considered in more detail via the 
coherence function, which we also use for t/t' extrapola- 
tion. We define this function as follows in order to remove 
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FIG. 5: The decoherence functions D{ky/-Fr, t/t') in the matter era from 1024'' simulations with s = 0, averaged over 3 
realizations. The raw data is highlighted by the lighter central region, with extrapolations to more extreme t/t' values and to 
lower ksfrr' are indicated by the darker regions. 



the equal-time fcr dependence from the unequal-time re- 
sults: 



b{k\RT' ,tIt') = 



C{kV^,T/T') 



C{kT,l)C(kT',l) 



(28) 



where the modulus in the square-root is relevant only 
for (7^2- This has the attractive feature of being equal 



to ±1 at equal-times and is positive at such times for 
all auto-correlations. Our results are shown in Fig. [SI 
Fig. |6l and also in Fig. |9l The small-scale behaviour 
can be understood by considering that the network can 
quickly decohere when coarse-grained on sub-horizon 
scales 1/fc ^ ^ since the relativistic strings must simply 
travel ^ 27r/fc. As a result the coherence function decays 
by |fcr — /cr'l ~ 27r. Figure [6] highlights this decay, and 
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FIG. 6: Slices through the decoherence functions D at small 
scales for the matter era. In order to demonstrate the ap- 
proximately constant form at these scales when plotted as a 
function of A;(r — r'), results are shown for the three values 
kij -\- r')/2 indicated by the legend. The latter quantity is 
perpendicular to k(T — t') in the kr — kr' plane (see also 
Fig. [9| and simplifies to fer for the equal-time case. 



that the form of I? as a function of A:(r — r') is scale- 
independent on small scales. On the other hand, when 
plotted as a function of t/t' , as in Fig. [51 the equal-time 
ridge becomes increasingly sharp on these scales. 

On super-horizon scales the decay is much slower, tak- 
ing of order a Hubble time: t/t' ^ 2. For these longer 
wave modes the coarse-grained regions each contain a 
number of horizon volumes. A particular horizon vol- 
ume that is initially over-dense will become under-dense 
due to the stochastic string dynamics in a time that is of 
the same order its initial size, i.e. '^ t' . Therefore the 
averaged properties of the large coarse-grained region de- 
cohere on this time-scale and the decoherence functions 
for large scales show decay at a fixed value of t/t', inde- 
pendent of scale. 

We extrapolate our UETC scaling function results to 
greater time ratios in the auto-correlation cases by first 
noting that the form of the coherence function at fixed 
k\/ tt' is approximately Gaussian in log(T/r'), but with 
the width dependent on kyprr' . For super- and near- 
horizon scales we hence take each kyrr' value in turn 
and match a Gaussian profile exp[— log {t/t') /2(t'^] to 



the data at the most extreme t/t' value and then use 
that profile to extrapolate to larger time ratios. On the 
other hand, for sub-horizon scales the decay is very rapid 
and we simply extrapolate with zeros. The results of 
these extrapolations can be seen in Figs. |4]and[5] 

For the cross-correlator C'j'2 on the other hand, our 
results show that the r/r'-profile on horizon and super- 
horizon scales is complex and non-Gaussian, making reli- 
able extrapolation impossible. For this correlator we use 
the Gaussian fits to Cfi and C22, with widths an and 
(T22, to provide information on the likely decay timescale. 
We then fit a Gaussian of width (J12 = {(Ju + oi-i)/2 but 
with free mean and normalization to the Cy^ data from 
the most extreme time ratios, and then employ that to 
yield extrapolation. 



D. Extrapolation of UETCs to sub-string scales 

Figure [7] shows a plot of kTC\^ for t' = t from simula- 
tions with s — Q. From the upper pane it can be seen that 
this quantity is approximately independent of fcr on sub- 
horizon scales until /cr ~ 200 at the earliest time shown 
or until about kT ^ 400 at the latest time plotted. The 
plateau highlights an important point: that Cfi drops 
off as « 1/kT on sub-horizon scales, which matches our 
basic expectations as well as similar measures from NG 
simulations [64]. However, on smaller scales there is a 
sudden attenuation of power, and at a fcr value that is 
increasing with time. The lower pane clarifies the later by 
re-plotting this against fc rather than fcr, showing that 
this is occurring at a fixed comoving scale of fc sa 20o. 
This corresponds to a few times the string width, which 
is a fixed comoving width in this s = simulation. 

We find evidence that all equal-time scaling functions 
vary as approximately 1/ fcr on small scales (except in the 
vector case for which [kT^C^ varies roughly as 1/fcr). 
Figure [5] shows power-law fits to fcrC (and {kT)^C in the 
vector case) over the range fcr = 30 — >■ 100, which lies 
between the interesting effects on horizon-scales and the 
artifacts of the string width on very small-scales. While it 
would be desirable to have a further order of magnitude 
via which to confirm the behavior in this regime, it is 
clear in all cases that a substantial improvement in our 
estimate of the scaling functions on small scales would 
be arrived at by extrapolating this trend down to scales 
near and below the string width. 

We extrapolate the unequal time scaling functions to 
small scales using our knowledge of the attenuation expe- 
rienced at equal-times for both r and r'. First we define 
the attenuation level R as the ratio of the extrapolated 
equal-time scaling function and the measured version at 
time r: 



R{k,T) 



Q{kT) 



a(fcT,l) 



(29) 



where we have added a subscript to C to indicate that 
the measured correlator does not scale exactly and so is 
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FIG. 7: Results for krCfi from s = simulations, plotted 
against kr (upper pane) and k (lower pane). The 5 lines cor- 
respond to equally spaced times between r = 150 and 3004>q^ . 
In the upper pane the 5 lines are indistinguishable at low kr 
due to the observed scaling, but scaling is broken on small 
scales where lines progressively move to the right. In the 
lower pane the lines are indistinguishable at high k, showing 
that on such scales krCfi is then simply a function of k, ie. 
the comoving length scale, and indicating that the attenua- 
tion of the signal occurs on scales a few times the comoving 
string width, which is constant in this s = case. 
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FIG. 8: Results for all five equal-time scaling functions mul- 
tiplied by kr (or (kr) in the vector case), with power-law 
fits over the range fcr = 50 — ^ 100 (between the two vertical 
lines) . The uppermost line is for Cfi , then of the middle pair 
of lines the black one is C22 and the gray one is 1(5121, and 
finally of the lower pair the black one is (7v and the gray one 
if (TT. A horizontal fit line would indicate that C (or {krYC^ 
in the vector case) varies as 1/fcr. 



not only dependent upon fcr, while p and Q are constants 
found from the power-law fit. We then extrapolate the 
UETC scaling function as: 



7t') == y/R{k,T)R{k,T') Cr- (30) 

This is based on the correlators being quadratic quanti- 
ties and includes an appropriate compensation factor for 
each of the two times involved, which in practice must 
be nearly equal for the scaling function to be signifi- 
cant, while in the equal-time case this simply returns 
the power-law fit. 



E. CMB calculations and UETC eigenfunction 
decomposition for small scales 

As explained in Sec. HH h aving determined the UETC 
scaling functions C{k\/ tt' ^ t Jt'\ we decompose them 
into a sum of terms involving functions of a single variable 
Cnikr), which then act as the sources of metric pertur- 
bations in the CMB calculations. Each function hence 
corresponds to one term in the sum for the CMB tem- 
perature or polarization power spectra. 

The following discussion is simplified if we discuss ini- 
tially only the tensor case, for which we require a single 
UETC that is symmetric under the exchange of r and 
r'. Since we can only represent Cnikr) numerically at 
discrete fcr values Ki (with i = 1, . . . , M), then this step 
boils down to the decomposition of a real and symmetric 
M X M matrix Ctj = C{y/KiKj, Hi/Kj), the eigenvectors 
of which form an orthonormal set. It can hence be seen 
that the decomposition: 



M 



Gi 



Ea. 



(31) 



is equivalent to determining the eigenvalues A„ and eigen- 
vectors Cni of the matrix C'ij. In the scalar case, the 
situation is more complex but the discussion proceeds 
similarly"^, while in the vector case we apply the decom- 
position to k^^/rr'C. 

In BHKU we found that accurate results over the rel- 
evant angular scales were obtained with M = 512 and 
Km = 200, while only the terms with the 128 highest 
|A„| values were needed. However, here we wish to in- 
clude the high-fcr tails of the UETCs out to fcr - 5000, 
but the narrow width of the equal-time ridge requires 
An « 1 and hence we would require M ~ 5000. Addi- 
tionally, the small amplitude of the tails implies that their 



^ In the scalar case there are two source functions S$ and 5* 
between which there is a finite correlation. To deal with this 
we form a 2M X 2M symmetric matrix by tiling Cfi, Cf2j f^fl 
and C22 ^iid then take S$ to be the first M elements of the 
eigenvector and S-q/ to be the second half (see BHKU). However, 
this barely changes the present discussion of small scales 
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signal is likely to be contained in eigenvectors with very 
low eigenvalues, and therefore we would need to include 
all terms in our CMB calculations. Due to the nature 
of the sources, each Einstein-Boltzmann integration is 
much slower than the corresponding primordially-seeded 
calculation required for CMB predictions from inflation 
and hence this process would be particularly time con- 
suming. Further, the contributions to the CMB power 
spectrum in our target range 2 < ^ < 4000 from ex- 
tremely high kr are minor, while our knowledge of the 
UETCs on such scales is only via the above power-law 
extrapolations. 

We proceed instead by performing the decomposition 
in such way that each "eigenvector" is localized in kr, 
since we then have an immediate understanding of how 
it contributes to the CMB power spectra. Firstly we ar- 
range for all of the dominant horizon and super-horizon 
power [kr < 100) to be contained within a particu- 
lar set of eigenvectors, which then completely dominate 
the CMB temperature power spectrum for multipoles 
i < 1000, but their contributions decay for smaller angu- 
lar scales. Secondly, the simple form of the UETC scal- 
ing functions on sub-horizon scales allows us to obtain 
a knowledge of the CMB contributions from extremely 
high kr values by combining our calculations for moder- 
ate kr values with approximate scaling laws, as explained 
momentarily. 

The fcr localization is achieved by noting that the rapid 
decay of the scaling functions for unequal times means 
that Cij can be approximated as the sum: 






(32) 



where, as indicated in Fig. [HI CJ" is finite only for 
<^m < hi < 5m+2^ with the 5n chosen to yield an ar- 
ray of overlapping matrix blocks which cover the im- 
portant regions of the kr-kr' plane: horizon scales and 
the sub-horizon equal-time ridge. In the regions of over- 
lap, the UETC power is shared between the two ma- 
trix blocks such that the contribution from a given block 
varies smoothly from zero at the extremes of the corre- 
sponding kr range, up to full at the block centre. As a 
result of this construction, each component matrix CJ" 
has eigenvectors c™" which have the desired locality in 
kr, being finite only for (5m < * < 3„i+2, while there are 
Sm+2 — Sm eigenvectors for each block. We chose the 
Sm values and kr spacing such that the first block is fi- 
nite only for kr < 100, while the subsequent blocks have 
width A(fcT) = 20, as indicated in Fig.[S] 

Importantly, the content in all of the higher blocks is 
of the same form: a ridge of given width with a cen- 
tral height that decays by a fraction ~ Afcr/fcr across 
the block (in addition to the decays required to share 
power between overlapping sub-matrices). So long as 
Afcr << kr, then the ~ l/fcr power-law decay has mini- 
mal impact and the differences between the higher blocks 
are effectively just translation in kr and a change in nor- 
malization. That is, the eigenvectors are effectively trans- 
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FIG. 9: The block-wise coverage of the kr — kr' plane. A large 
primary block contains the dominant UETC data at near- and 
super-horizon scales (kr < 100). Then smaller blocks of side 
A(fcr) — 20 cover the diagonal, with overlapping to aid cover- 
age. The gray-scale image shows the coherence function Dfi 
to illustrate that these patches cover the region over which 
the UETC is significant, with light colors indicating higher 
values. The dark outer region indicates the region for which 
we extrapolate the coherence function, as in Fig. (5] 



lated in kr while the eigenvalues absorb the normaliza- 
tion change. 

When our modified version of CMBEASY is applied to 
the eigenvalues and eigenvectors from the higher blocks, 
the temperature and polarization power spectra returned 
also have effectively the same form, except for the ap- 
proximate re-scaling: 



PI' 



£{£ + 1)C^ 



fie/^im) 



J 



(33) 



where k,„ is the central fcr value of the mth block while 
/? w 2 is a constant to be determined. This is illus- 
trated for the vector mode in Fig. [101 It can also be seen 
that the contribution from a given block rises as approx- 
imately oc € up to a plateau for £ w Km — > 200^™ before 
rapidly decaying for higher £. The key point is that given 
this scaling, we can find an approximation to the CMB 
power spectrum contribution from fcr values beyond the 
maximum value for which we perform full calculations 
Kfuib by simply summing over many such re-scaled forms 
until the desired accuracy is reached. The total power 
spectrum is hence: 



T', = g r^- 



E 



f{i/Km) 



(34) 
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FIG. 10: Contributions to the vector component of the tem- 
perature power spectrum from discrete kr ranges. The high- 
est hne shows the dominant contribution from horizon and 
super-horizon scales. The contribution from first sub-horizon 
matrix block is shown in the middle of the plot, while the ap- 
proximate scaling law is indicated by the high kr lines, which 
are lowest in the figure. 



FIG. 11: The build up of the high I scalar mode tempera- 
ture power spectrum as the maximum included fcr value is 
increased. From bottom to top in the plot these correspond 
to the fcr < 150, fcr < 300, fcr < 550 and finally fcr < 5000. 
Results shown by thin solid grey lines are from full CMB cal- 
culations, while the thick grey dashed and thick black solid 
lines correspond to results obtained using the rescaling prop- 
erty for fcr > 150. The power-law extrapolation is shown by 
the thin dashed line and is a reasonably accurate description 
of the final kr < 5000 result for £ > 3000. 



which is shown in Fig. [TT] for the scalar mode. Further- 
more, since Vp decays rapidly for £ > 200Km, then we 
may extrapolate the sum for £ — > oo with mmax — >■ oo by 
approximating it as an integral, which then varies as: 



IV. CMB RESULTS 
A. Temperature Power Spectrum 

1. Results 



ex 



£0- 



(35) 



independent of the detailed form of the function f{£). 
For /3 sa 2, this yields roughly the l/£ form expected at 
very small angular scales due to the GKS effect ,33], as 
will be discussed in Sec. IIVA2I For the scalar, vector 
and tensor modes respectively we find P — 2.5, 1.9 and 
1.7 and we hence tentatively predict that the CMB con- 
tributions from each vary as £~^-^, £~'^-^, £~^-'^ at very 
high £ (see next section). 

Note that the use of this approximate re-scaling prop- 
erty is applied only to fcr values for which our knowledge 
of the UETC scaling functions is arrived at via power- 
law extrapolation. That is, we perform full CMB cal- 
culations for the range of kr for which the simulations 
provide direct information, and it is fortunately the case 
that these provide the overwhelmingly dominant contri- 
bution for scales £ < 1500, while in the present article we 
are limiting ourselves to .^ < 4000 and hence the extrap- 
olations made are fairly reliable. 



We present our final CMB power spectrum results for 
cosmological parameters: h — 0.72, Q,bh? — 0.214 and 
r^A = 0.75, which match those used in BHKU and the 
central values from non-CMB measurements [65 - 67] . We 
additionally assume an optical depth to last-scattering 
of ?7 = 0.1, which again matches BHKU. For the case 
of the temperature power spectrum these are shown in 
Fig. [T^ The form is essentially that found in BHKU: at 
low £ there is a roughly ex log£ rise up to a broad peak 
between .£ « 30 -> 700 (75% of peak), with the peak itself 
at £ « 400. At greater £ the spectrum decays, initially 
as roughly l/i^, but this then slows to roughly \/£ for 
£ > 3000, which is more clearly evident in Fig. [T3l 

The reason for this change in behaviour can be seen 
in Fig. I14[ in which we compare our results to those in 
which the cosmic string sources are artificially zeroed for 
times prior to recombination. This zeroing can be seen to 
remove from the signal a contribution that dominates for 
£ « 500 — > 2000 over a smaller signal that is sourced at 
later times and that has the approximate 1/ £ form. That 
is, the early-time signal is Silk-damped at small-scales 
such that the late-time signal is revealed and dominates 
for £ > 3000. Note that this decomposition into pre- and 
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FIG. 12: The CMB temperature power spectrum determined 
from s = simulations and incorporating estimated UETC 
power for fcr < 5000. The plot shows the total (thick line) 
plus the decomposition (thin lines) into scalar (S), vector (V) 
and tensor (T) modes. 
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FIG. 14: The string contribution to the temperature power 
spectrum when including the string sources only after recom- 
bination (thin) or at all times (thick), with the difference be- 
tween the two additionally shown (dashed). Recombination is 
not instantaneous and additionally the string source functions 
must be temporally differentiable and thus must be gradually 
switched on, hence this decomposition is ambiguous and there 
is an artificial reduction in the amplitude of the power-law at 
high H for the post-recombination results. 



is also interesting to note what expectations we have for 
I > 4000. As discussed in the previous section, the 
^ -> oo limit of the Eq. (IMl) sum suggests the contri- 
butions vary as i~P, where p — 1.5, 0.9 and 0.7 for 
the scalar, vector and tensor modes respectively. How- 
ever we caution against large extrapolations using these 
power-laws since our simulations do not have the dynam- 
ical range required to yield confident knowledge of the 
UETC power-laws at high fcr, particularly in the vector 
and tensor cases, and that uncertainty feeds through to 
the present power-laws (see Fig. [5]). 



FIG. 13: Results for l^'Ce, highlighting the form of the power 
spectrum at high i. The plot shows the total (thick line) 
plus the decomposition (thin lines) into scalar (S), vector (V) 
and tensor (T) modes. Additionally, the right of the figure 
shows the results of our tentative power-law extrapolations to 
I > 4000. 



post-recombination is ambiguous since decoupling is not 
instantaneous and additionally we require the scalar and 
vector sources to be temporally differentiable, implying 
that we must have a gradual transition from zero to the 
nominal source value. These are responsible for the small 
difference in amplitude between the two power-laws seen 
in the figure. 

Although the aim of our calculations was merely to 
obtain results at ^ < 4000 that are accurate enough to 
allow reliable comparison against observational data, it 



2. Comparison with results from Nambu-Goto strings 

Firstly our results obtained for small angular scales 
are in broad agreement with analytical expectations from 
the GKS effect using a Gaussian model of a Nambu-Goto 
string network [33l] . which predicts a 1/i dependence in 
the small-scale limit. That model would have yielded 
ETC functions behaving as l/fcr at small scales, which 
is not precisely as seen in our simulations. 

Further, this approximate dependence has also been 
confirmed in work based upon high resolution Nambu- 
Goto simulations by Fraisse et al. (2008) [39,]. While the 
Nambu-Goto simulation result stems from a method that 
does not include recombination and hence only includes 
the effects of strings after last-scattering, their results re- 
veal a l/£^'^ variation in the power spectrum over the 
range 400 < ^ ^ 10^. Their calculations are hence most 
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comparable^ to our results in Fig. [Ml when we artificially 
zeroed our string sources until decoupling, which also in- 
dicate a power-law dependence for (. > 400. Importantly 
our full results demonstrate for the first time at which 
angular scales this « l/£ dependence is valid, namely 
£ > 3000 — scales finer than a few arc minutes. 



At angular scales near the peak of the spectrum {£ ~ 
400) there are no recently published results from Nambu- 
Goto simulations against which we can compare our re- 
sults, with Contaldi et al. (1999) [40| providing the 
only example of such work. However, this employed 
Minkowski space-time simulations plus now out-dated 
cosmological parameters and therefore a detailed com- 
parison is not appropriate. The basic form of our spec- 
trum is in good agreement with their results, although 
their power spectrum shows a shift to higher £ compared 
to those presented here. Further, their results suggest 
that Nambu-Goto networks yield CMB predictions with a 
larger overall normalization. This is also seen by compar- 
ing to power spectra calculated from FRW Nambu-Goto 
simulations and modern cosmological parameters, but 
which are valid only for discrete multipole ranges (and 
miss the i ^ 400 peak) [H HI]. Nambu-Goto strings 
would require G/x w0.7xl0~^in order to fit observations 
at £ < 10 [H, while we would require G^i = 1.8 x 10"^ 
to match the WMAP5 resuh at £ = 10. This is not 
unexpected given the simulation results and, as shown 
in Table U Nambu-Goto calculations yield higher string 
densities than field theoretic ones, raising the power spec- 
trum normalization and shifting power to smaller scales. 
We believe this density difference is because field theory 
simulations provide a decay channel, namely radiative 
decay, which is not included the Nambu-Goto codes. De- 
spite this, the Fraisse et al. result for £{£ + l)Cf/(G/i)^ 
at £ = 4000 is ~ 25, which compares favourably to our 
result of 20 at this multipole. This may be due to am- 
biguities associated with their non-inclusion of recombi- 
nation rather than a real agreement in the amplitude of 
the power-law. 

As discussed in the introduction, the USM is a com- 
putationally cheap means of estimating the CMB power 
spectrum, and this has also been used to study the string 
contribution at sub-WMAP angular scales with the USM 
parameters set to ap pro ximate a Nambu-Goto network. 
The published work [43| highlighted a l/£'^ dependence 
near £ ~ 2000, but a more recent look at greater £ 
prompted by our results does indeed reveal an approxi- 
mately l/£ variation in the USM power spectrum above 
£~3000 [69|]. 




10 100 



1000 



2000 



3000 



4000 



FIG. 15: A comparison of three contributions to the CMB 
power spectrum: the inflationary scalar contribution (grey, 
solid), the string contribution (black) and the SZ effect (grey, 
dashed). The inflationary scalar mode and /lo « 1 string 
contribution are both normalized to match the WMAP data 
jj,] at £ = 10, while the /lo ~ 0.1 string contribution is set 
to yield 10% of the observed power at this multipole (the 
approximation signs are because /lo is defined relative to the 
total theoretical power spectrum, not the observations). Two 
normalizations are shown for the SZ contribution, one for each 
of the QUAD bands: 100 and 150 GHz. 



3. Comparison with the contribution from inflation 

It is of course useful to compare the cosmic string 
power spectrum with that from infiation, to which it 
should be added^. This comparison is eased if we give 
the string contribution an artificially high normalization 
such that it is equal to the infiationary contribution at 
i = 10 and then set both to be equal to the value observed 
by WMAP at that multipole, as in Fig.[T5l It can be seen 
that while the string contribution approximately tracks 
the troughs in the infiationary contribution for £ < 1500, 
the exponential suppression of the infiationary contribu- 
tion at high-^ means that the string component grows 
in relative size, and dominates for £ > 2000. Switching 
to a more realistic string contribution, with one-tenth of 
the previous normalization, this domination is delayed 
until £ « 3500, but importantly the fraction of the to- 
tal theoretical spectrum due to strings fi increases from 
/i500 ~ 0.1 to /3500 = 0.5 and therefore accurate data at 
£ > 2000 should in principle be highly sensitive to cosmic 
strings. 

However, these are the same angular scales for which 



"* Note that in contrast with Fraisse et al., even these results still 
include the effects of perturbations induced in the matter by the 
string sources and do not make use of the flat-sky approximation, 
hence they are still not including exactly the same CMB physics. 



^ We calculate the inflationary contribution using CAMB [7Cl |. 
with the same parameters as used in the string calculations and 
also rig = 1.0, while assuming negligible contributions from pri- 
mordial gravity waves. 
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the Sunyaev-Zel'dovich (SZ) effect [71| begins to make a 
significant contribution to tlie temperature power spec- 
trum at certain observational frequencies. Tlris contribu- 
tion results from the distortion of the black-body spectral 
profile as the CMB photons pass through galaxy clusters 
and Compton-scatter off hot electrons. In Fig. \TE \ w e 
plot the predictions made by Komatsu and Seljak [72l |. 
normalized^ for the two frequency bands observed by the 
QUAD project ^^, namely 100 and 150 GHz. It can 
be seen that the string contribution with /lo ~ 0.1 is 
likely to be shrouded, even at high £, when measurements 
are made at low frequencies, but in observational bands 
near 220 GHz, the SZ effect is suppressed. Unfortunately 
higher frequencies have a greater contribution from un- 
resolved point sources J30!] and in practice observations 
at a number of frequencies will be required in order to 
understand the frequency dependent contributions and 
remove them to yield a high sensitivity to the frequency- 
independent cosmic string component from temperature 
power spectrum measurements alone. 



B. Polarization Power Spectra 

CMB polarization anisotropics differ from those for 
temperature in that they are created almost exclusively 
at recombination. Since the polarization is caused by 
Thompson scattering in the presence of a quadrupole 
anisotropy, the universe must be ionized in order to cre- 
ate it, but with Thompson scattering still weak enough 
for it not to have suppressed the quadrupole too heavily. 
Polarization is hence created only very close to recom- 
bination, and also weakly during the recent re-ionized 
epoch. As a result, the polarization signal comes primar- 
ily from horizon-scale UETC power impacting at times 
close to recombination. The contributions from high fcr 
near last scattering are outside our window of interest 
£ < 4000 and are suppressed by Silk-damping. Further, 
the contributions from high kr at recent times are neg- 
ligible because there is little UETC power at such scales 
and the re-ionization contribution is secondary in impor- 
tance. Our present results, shown in Fig. 1161 therefore 
add little new information for CM B p olarization from 
strings beyond our previous results [22, l2J] . 

However, for completion, it should be noted that while 
the string contribution is sub-dominant for all scales in 
the E-mode polarization spectrum (EE) and the cross 
correlation of the E-mode with temperature spectrum 
(TE), it may dominate the B-mode polarization spec- 
trum (BB) for scales i < 400. This is possible because 
the inflationary scalar mode contributes to BB only via 
weak gravitational lensing, which converts EE power into 
BB. Furthermore, our results indicate that the string con- 






o 

+ 




The plotted SZ spectrum was calculated for fli,h^ = 0.023 and 
(Tg = 0.8 and scales in normalization as roughly (Q^^h)^a'l [7211 . 



FIG. 16: A comparison between the string contribution 
(thick) to the CMB temperature and polarization power 
spectra with the adiabatic scalar contribution from inflation 
(thin). The normalization of the string component is set at 
/lo « 0.1, while the inflationary scalar contribution uses the 
same settings as in Fig. [TS] For the BB spectrum we also plot 
the possible contribution from the inflationary tensor mode 
(dot-dashed), normalized at its 95% upper limit |73j ]. 



tribution peaks at significantly lower £ than the weak 
lensing contribution, which should readily prevent con- 
fusion between the two, while the weak-lensing contam- 
inant can also be partially removed [TJ. Note that we 
find a BB peak from strings that is at larger scales than 
is seen in USM results when that model is set to approx- 
imate Nambu-Goto simulations [4^, as expected given 
our lower string densities, and that the USM results sug- 
gest more confusion with the weak-lensing signal. It is 
also possible that the inflationary tensor mode may make 
a sizable contribution to the BB spectrum at very large 
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angular scales, as shown in Fig. [TBI but its normalization 
depends on the details of the inflationary model and is 
poorly constrained by current data. A wealth of B-mode 
data will soon be available and promises to be highly 
sensitive to cosmic strings [2J, |4y, [7J| . 



V. CONCLUSIONS 

We have presented the first calculation of the cosmic 
string CMB temperature power spectrum contribution 
that is accurate over the multipole range i — 2 ^^ 4000, a 
range that encloses the scales probed by the Planck satel- 
lite (^ = 2 — > 2500) and additionally scales at which sub- 
orbital data are becoming increasingly accurate. Further, 
we show that aX i ^ 3000 there is a knee in the power 
spectrum result due to the exponential decay of the early 
time contribution at small-scales, which then reveals the 
post-recombination component. This late-time contribu- 
tion varies as roughly l/£, which is the basic expecta- 
tion for the GKS effect [33!|. Our results yield values 
for the power law exponent for each of the scalar, vector 
and tensor contributions to the overall high-£ behaviour 
and can in principle be used to estimate the temperature 
contribution from strings for £ > 4000. However this 
extrapolation must be performed with caution since our 
simulations do not have the dynamic range to yield these 
power laws with great confidence. 

Our results indicate that the size of the cosmic string 
signal in the temperature power spectrum, relative the 
adiabatic inflationary component, is roughly 10 times 
greater at £ w 3500 than it is at ^ w 10 or £ w 400 - 1500 
and therefore small scales are particularly important 
for cosmic strings. While it is true that other con- 
tributions also become significant at small scales, for 
example the Sunyaev-Zel'dovich effect, such effects are 
frequency-dependent and hence can be identified and 
subtracted. Further, our results are not limited to tem- 
perature anisotropies but include polarization also, which 
is of great future importance for strings in the case of the 
B-mode. We present a comparison of our temperature 
and polarization results to the latest data in a separate, 
shorter article [73]. 

Our simulations also yield greater knowledge of the 
scaling properties of Abelian Higgs string networks, with 
improvements in the accuracy of results such as the scal- 
ing density and the UETCs. Additionally, we have mea- 



sured the coherence function for Abelian Higgs strings 
for the first time, the form of which is important for cal- 
culations of this kind. 

That our results stem from the Abelian Higgs model 
means, of course, that they are not necessarily accu- 
rate for cosmic superstrings (see eg. [ifll for a review). 
These superstrings may have intercommutation probabil- 
ities significantly lower than that for gauge strings and 
additionally may form Y-shaped junctions, which do not 
form in the Abelian Higgs model for the parameters cho- 
sen here. While small-scale structure on the strings has 
been shown to lessen the impact of the intercommutation 
probability ^^, the effect of Y-j unctions on the network 
properties is highly uncertain. Both effects are likely to 
increase the string density and decrease the inter-string 
distance, resulting in a shift of the peak in the CMB 
signal to greater £ values. This may therefore enable su- 
perstrings to be distinguished from conventional cosmic 
strings should a string component in the CMB be de- 
tected in future data. 

Finally, we note that the power spectrum is not a com- 
plete statistical description of the anisotropies that would 
be seeded by cosmic strings, since there would be a sig- 
nificant non-Gaussian character created by the GKS ef- 
fect. Higher order moments such as the bispectrum and 
trispectrum have been calculated |77H79{ , in addition to 
realizations of CMB maps [4J, [gg ■ While these calcula- 
tions are challenging and either do not include recombi- 
nation or are valid only for very limited multipole ranges, 
non-Gaussianity is an exciting channel for future cosmic 
string constraint or detection. 
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